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1. Introduction 

Over the past decade, MINOS has been successfully measuring fundamental neutrino oscilla¬ 
tions parameters. Using the NuMI beam ||T]] and its two-detector design, a Near Detector (ND) at 
the beam source at Fermilab and a Far Detector (FD) 735 km away in the Soudan Mine in Min¬ 
nesota [Q], MINOS was able to cancel out many systematic uncertainties, for example associated 
with Neutrino Flux [|^, and push the boundaries of the measurements that could be made |Q]-[^]. 

MINOS took data between May 2005 and April 2012, mostly in the so-called Low Energy 
(LE) beam configuration, with the neutrino spectrum peaked at approximately 3 GeV, close to the 
standard 3-flavour oscillations maximum at about 1.6 GeV in the MINOS Ear detector. This al¬ 
lowed MINOS to carry out precise measurements of the muon neutrino disappearance oscillations 
parameters. MINOS was then reborn as the ongoing MINOS-i- experiment |^, and is now taking 
data in the NOvA-era NuMI beam, which started producing neutrinos in September 2013 and is 
running at a higher energy, optimal for the NOvA experiment ||^. This so-called Medium Energy 
(ME) beam is peaked at approximately 7 GeV, relatively far from the standard 3-flavour oscillations 
maximum, however, on the other hand, MINOS-i- is currently the only neutrino experiment in the 
world seeing a wide-hand on-axis beam, with a wider neutrino energy range than other accelerator 
neutrino experiments such as NOvA and T2K, and with significantly more events per unit of expo¬ 
sure due to being on-axis. MINOS-i- has therefore a unique potential to look for new physics and 
exotic phenomena such as sterile neutrinos, Non-Standard Interactions (NSI) 0, and large extra 
dimensions. This proceedings paper summarises some MINOS and MINOS-i- results as of October 
2014. 

1.1 The MINOS(-i-) Experiment 

The MINOS(-i-) experiment consists of two detectors, a Near and a Ear Detector. The detectors 
are in the path of the NuMI neutrino beam. The ND is located at Eermilab, 1.04 km from the NuMI 
beam source, and its purpose is to measure the neutrino beam before oscillations have occured. The 
ED is located ~700 m underground in the Soudan Mine, and measures the neutrino beam after the 
neutrinos have oscillated. 

In order to generate neutrinos, the NuMI beam sends bunches of 120 GeV protons onto aim 
long graphite target. This generates hadrons which are focused by two magnetic horns and decay 
into secondaries and neutrinos in a 675 m long decay pipe. The secondaries are then removed from 
the beam by a hadron absorber at the end of the decay pipe and a “muon shield” consisting of 240 m 
of dolomite rock on the way to the MINOS ND. This leaves a neutrino beam that consists mostly 
of muon neutrinos (91.7%), with an admixture of muon antineutrinos (7%), and a small electron 
neutrino contamination (1.3%). The NuMI beam can be modified to produce an antineutrino en¬ 
hanced beam (rather than the standard neutrino-dominated beam) by reversing the horn current and 
thus focusing negative hadrons. 

Both MINOS(-i-) detectors are large magnetized iron scintillator tracking calorimeters opti¬ 
mised to measure muon tracks from muon neutrino charged current (CC) interactions and to sep¬ 
arate muon neutrinos from antineutrinos based on track curvature. They are built out of alternate 
planes made of steel and scintillator. The scintillator planes are attached to the steel planes and 
consist of long scintillator strips with embeded wavelength shifting fibers that are read out by pho- 
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Figure 1: Pictures of the two MINOS+ Detectors. The Near Detector is smaller, at about 1 kT mass and 
15 m length. The Far Detector is very large at 5.4 kTon mass and 30 m length, and has a veto shield above 
it to identify cosmic ray muons. 


tomultiplier tubes. Subsequent scintillator planes are positioned orthogonally to each other for 3-D 
reconstruction of events. The steel planes serve as a large detector mass to increase the likelihood 
of neutrinos interacting in the detector. Figure [T] shows pictures of the MINOS Near and Far de¬ 
tectors. The ND is 15 m long, 3.8 m tall and 4.8 m wide. It weighs approximately 1 kT and has a 
“calorimeter” region where each steel plane is instrumented with scintillator, and a “spectrometer” 
region where only every fifth plane is instrumented; there are 282 steel planes, whether instru¬ 
mented or not, in total. The FD is 30 m long, 8 m wide and 8 m tall. It weighs 5.4 kT and consists 
of 486 fully instrumented steeFscintillator planes. It also has a veto shield to identify cosmic ray 
muons. 


1.2 Neutrino Oscillations 


There are three known active neutrino flavours, the electron, muon, and tau neutrino {Ve,VpL,Vz)- 
Each neutrino also has a corresponding antineutrino. Neutrinos interact only via the weak force 
(and gravity), and therefore very rarely, and have been observed to oscillate between the flavours 
as they travel. Neutrino oscillations can be described by Equation |1.1| : 


VVr 


= u 



(1.1) 


where U is the Pontecorvo-Maki-Nakagava-Sakata (PMNS) matrix [[T^]-[[T^ which has a 
parametrisation of the form: 


c/ = 


C12C13 S12C13 il3e“ 

—il2C23 — Ci2i23-«13e'^ C 12 C 23 — ^ 12^23^13^'^ i23Cl3 

■Sl2i23 — Cl2C23-Sl3e'® ^Cl2i23 " ■«12C23‘S13'^'^ C 23 C 13 


( 1 . 2 ) 


In the above, are the neutrino flavour eigenstates, Vi, 2,3 are the neutrino mass eigenstates, 
Cij = cos Qij, Sij = sin djj, with dij being the mixing angle, and 5 is the charge-parity (CP) violating 
phase. 

Eor a two-neutrino approximation where the difference in mass between two of the mass states 
is very small, the probability that a muon neutrino will have oscillated to a tau neutrino as it prop- 


3 










MINOS and MINOS+ Results 


Anna Holin 


agates is to first order given by: 


Pv.,^vr = sin^20sin^ 

M ^ 


l.llNirPL 


(1.3) 


where L is the baseline in km, E is the neutrino energy in GeV, and is the mass squared 
difference between the two effective flavour states in eV'^ jc^. Three-flavour neutrino oscillations 
are actually parametrised using two mass squared differences - (solar sector) and Ani 32 (at¬ 
mospheric sector), three mixing angles - dn, 0 i 3 , and 623 , and a Charge Parity violating phase - 
5 . As the 013 angle has now been measured by the Daya Bay [|T^], Double Chooz |[I^ and RENO 
[]_5] experiments to be ss 9°, it can no longer be treated as negligible, and the full three-flavour 
framework has been adopted by MINOS for oscillation measurements. Currently the outstanding 
unknowns in the standard neutrino oscillations framework are the neutrino mass hierarchy (whether 
the mass state m 3 is heavier or lighter than the m\-m 2 mass state “doublet”), the CP violating phase 
8 , and the precise value of the mixing angle 623 , including the octant of the sin^ B 23 term. 

1.3 Sterile Neutrinos 

Though LEP has measured the number of light neutrinos to be 2.984±0.008 [ 16] ^ there is still 
the possibility that there are sterile neutrinos that do not interact via the weak force. There could 
be a fourth (or more) neutrino mass state that oscillates together with the standard active neutrinos 
but cannot be detected via normal charged current (CC) and neutral current (NC) interactions. Eor 
example, assuming a simple 3-1-1 model, this could lead to the observation of additional neutrino 
disappearance depending on the value of the third mass squared difference term Am 43 . There would 
also be additional mixing angles 614 , 624 and 034 . 

Sterile neutrinos could explain some of the anomalies seen in the neutrino sector such as the 
ESND [p^ and MiniBooNE [fi^] results, but so far the evidence has been inconclusive with tension 


between appearance and disappearance results. 


2. Results 


2.1 Standard Oscillations Three-Flavour Results 

New results are presented for the combined three-flavour analysis using disappearance data, 
Ve appearance data, and an increased atmospheric event sample of 48.7 kT-years. As in the previ¬ 
ous results from this analysis ||^, the full MINOS beam data set is used with a total exposure of 
10.71 X10^° protons-on-target (POT) of neutrino-dominated mode beam data, and 3.36 x lO^*’ POT 
of antineutrino-enhanced beam mode data. The larger atmospheric event sample used here corre¬ 
sponds to a 28% increase on the previous analysis, which used 37.88 kT-years of data. 

Eor the beam data, CC interactions are selected by requiring that there be a muon track 


and applying a k-Nearest-Neighbor algorithm [ ]20| ] using various muon track variables related to 
energy deposition and topology of the track. Both a contained-vertex (fiducial volume) sample is 
selected for the beam muon neutrino events and a non-fiducial muon sample, with vertices outside 
the detector. A special selection is also applied to select a sample of muon antineutrino events 


'This number has since been revised slightly downward to 2.92±0.05 0 - 


4 









MINOS and MINOS+ Results 


Anna Holin 


Neutrino beam 
(10.71x10^° POT) 


MINOS PRELIMINARY 


Antineutrino beam 
(3.36x10^° POT) 



Contained-vertex 

MINOS data 
Best fit 

No oscillations 
□ NC background. 
Cosmic-ray 


Reconstructed (GeV) 


- 800 

- > 

_ 0 600 

-V 

Non-fiducial muon8- 

fi 

> 30 
0) 

G 


7 Contained-vertex 

80 

G60 

V 

Contained-vertex . 

i 

ri 

- §400 
0) 



s 



§40 

0) 


-iV 

: iIJ 200 



iS 10 



iS 20 

T 







jj 

L _ 


6 9 15 

Reconstructed (GeV) 


10 15 20 25 

Reconstructed E^ (GeV) 


8 10 12 14 
Reconstructed E^ (GeV) 


Atmospheric neutrinos {48.67 kton-years) 

-Contained-vertex muons_^ _ 


- Non-fiducial muons - 


60-V^ E. = i-3Qiy 

Events 

V E^ = 3-10 GeV 

15 

|io 

V E, = 10-30 GeV 
|i 

30 

H20 

0 

m 

10 

V P^i<10GeV 

1 . 1 

Events 

.V P. >10GeV 

y -0.5 A fi.S 

008(0^) 


-0.5 A 0.5 

003 ( 62 ) 


008 ( 62 ) 


' ' ' '-0'.5' ' ' ' A 

008 ( 62 ) 


1 ' ' ' '-0'.5' ' ' ' A ' 

008 ( 62 ) 

25. y E, = 1-3QeV 

20 

^ E^= 3-10 GeV 

8 

n E, = 10-30 GeV 

2 d 
w id 

^ P^<10GeV 

^ i 

20 

• « 15 

^ P.>10GeV 

- 0.5 0 0.5 

H 

(B 

10 

-0.5 0 0.5 

1 4 
LU 

-0.5 0.5 

b 

Sid 

LU 

7 

-0.5 ( 

• I 10 

1 ^ 

Q 

- 0.5 0 


008 ( 0 ^) 


005 ( 0 ^) 


008 ( 0 ^) 


008 ( 0 ^) 


008 ( 0 ^) 


Figure 2: Main Samples used in Combined Analysis. MINOS beam data are shown in the top row, and 
the combined MINOS and MINOS+ atmospheric data for contained-vertex and anti-fiducial samples of 
neutrinos and antineutrinos are shown in the middle and bottom rows. The distributions of atmospheric 
neutrinos are plotted as a function of zenith angle and separated into bins of reconstructed neutrino energy. 
The observed data (black points) are compared with the prediction for no oscillations in gray and the best fit 
three-flavour oscillations in red. The cosmic ray background is shown in blue while the NC background is 
shown in gray. 


in the neutrino beam, and to select muon antineutrinos in the antineutrino-enhanced beam. For 
the Vg beam data, CC events are selected using a library-event-matching (LEM) algorithm, which 
compares shower events to large samples of tens of millions of NC background and CC-Vg signal 
events in order to determine whether a given event is CC-Vg-like [|1]. To include atmospheric neu¬ 
trinos and distinguish them from cosmic ray background events in the combined analysis, several 
different categories of events are analysed (see [ ^T| ] for detail). The two main sample categories 
are contained-vertex events and non-fiducial (upward going or horizontal muon track) events. The 
FD veto-shield helps in reducing the background cosmic ray muon events by checking for coinci¬ 
dences between hits in the shield and the selected tracks. The contained-vertex and non-fiducial 
afmospheric neufrino samples are furfher splif info muon neufrino and anfineufrino samples using 
frack curvafure, and info differenf energy slices. Finally, an afmospheric shower selecfion is also 
included in fhe combined analysis. A summary of fhe muon neufrino specfra used in fhis analysis 
is shown in Fig. ^ 

In fhe final fif, fhe Am 32 , 623 , and 5cp paramefers are leff unconsfrained. 0 i 3 is fif as a nuisance 
paramefer and consfrained by fhe reacfor resulfs (sin^ 613 = 0.0242 ± 0.0025). The solar mixing 
parameters are fixed fo Anj 2 i= 7 . 54 x 10^^ eV^ and sin^ 012=0.307. Systematic uncerfainfies are 
faken care of by fhe inclusion of 32 differenf nuisance paramefers in fhe fif. Due fo fhe separating 
of fhe neufrino and anfineufrino samples, fhe combined analysis is sensifive fo fhe neufrino mass 
hierarchy, fhe value of dcp, the octant of 623 , and the value of 0 i 3 . The calculations of oscillation 
probabilities in this analysis use the full three-flavour PMNS mafrix calculafions, wifh fhe inclusion 
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Data Set 

No Osc. Fred. 

With Osc. Pred. 

Observed Events 

from beam 

3201 

2496 

2579 

from beam 

363 

319 

312 

Non-fiducial from beam 

3197 

2807 

2911 

Atm. contained-vertex and 

1414 

1024 

1134 

Atm. non-fiducial and 

732 

575 

590 

Atm. showers 

932 

877 

899 


Table 1: Summary of Numbers of Events for the Combined Three-Flavour Analysis. 
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Figure 3: Contours and Profiles from Fit to 48.67 kT-years of atmospheric data combined with disappear¬ 
ance and appearance data from the MINOS beam. Left panels shows 68% and 90% confidence limits in 
(Am 32 , sin^ O^i) calculated for normal hierarchy (top) and inverted hierarchy (bottom). Right panels show 
log-likelihood profiles for each hierarchy plotted for A/rr^ (top right) and sin^ 023 (bottom right). The best fit 
is indicated by the star and occurs in the inverted hierarchy at Amjj = —2.37 x 10^^ eV^ and sin^ 023 = 0.43. 


of the MSW [|22|]-|23] effect. For the beam events, the ND is used to predict the FD event spectra. 


A summary of the numbers of events selected in the combined three-flavour analysis is pre¬ 
sented in Table |^. The fit obtains a best fit value for the atmospheric mass squared difference 
of Am 32 = 2.34lQQg x 10^^ eV^, and for the mixing angle parameter sin ^023 = 0-431 qq4 in the 
case of the normal mass hierarchy. The fit calculates values of Am 32 = 2.371qq7 x 10^^ eV^ and 
sin^ 023 = 0.43 ^qq5 in the case of the inverted mass hierarchy. The 90% confidence limits (C.L.) 
for the sin^ 023 parameter are 0.37 < sin^ 023 < 0.64 in the case of the normal mass hierarchy, and 
0.36 < sin^ 023 < 0.65 in the case of the inverted mass hierarchy. The best-fit contours are shown in 
Fig. m and show a slight preference for the inverted hierarchy and the lower octant of 023. A com¬ 
parison of the MINOS(-i-) results with the T2K experiment results [|^] are shown in Fig. The 
MINOS(-i-) results represent the most precise measurement of Am 32 to date. In the future, it should 
be possible to combine MINOS(-i-) and NOvA data together (when the latter becomes available) 
to achieve even more precise results and cut into the 023 octant and neutrino mass hierarchy phase 
space. 
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Figure 4: Comparison of contours from the MINOS(+) combined analysis in A77i32/sin^ 023 space with T2K 
experiment results . 

2.2 New MINOS+ Data 

The upgraded NuMI beam restarted on the fourth of September 2013 in ME mode and data 
from a total beam exposure of 1.68x10^° POT is used for the MINOS+ spectra shown in this 
paper. The new M1NOS+ data are presented in Fig. || together with the ratios to no oscillations. 
The oscillations prediction comes from the best-fit results presented in the previous section. The 
MlNOS-i- data are consistent with the MINOS oscillation results. The neutrino data sample contains 
1037 events with an oscillated (unoscillated) prediction of 1088 (1255) events. The antineutrino 
sample contains 48 events, with an oscillated (unoscillated) prediction of 47 (52) events. 

3. MINOS Sterile Neutrino Results 

MINOS uses the simple 3-1-1 model for its sterile neutrino results. Essentially it looks for 
distortions from a three-flavour formalism. There are thi'ee different regimes in which MINOS is 
sensitive to a fourth sterile neutrino. For small Am 43 <0.1 eV^, in the “slow” sterile oscillations 
regime, spectrum disortions are expected well above the standard oscillations maximum in the 
ED, with no effect in the ND. For medium values of 0.1 eV^ < Am 43 < 1 eV^, in the “interme¬ 
diate” regime, sterile oscillations at the ED become fast and average out and result in a counting 
experiment, still without effect in the ND. For large A 7 M 43 >1 eV^ however, in the “rapid” sterile 
oscillations regime, there are rapid oscillations at the ED which average out, and in addition there 
is an oscillations effect in the ND which affects the extrapolation to the ED. 

In order to account for ND distortions of the spectrum due to sterile neutrino oscillations, 
MINOS looks at the Far/Near ratio for both the CC and NC events and fits the oscillated F/N ra¬ 
tio directly to the F/N data ratio using a function. It fits for several pai'ameters: ^ 111 ^ 2 , 623 > 
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Figure 5: Reconstructed Energy Spectra and Ratios to No Oscillations for MINOS+ Far Detector Selected 
Events. Predictions are shown for no oscillations and for oscillations at the MINOS(+) best fit. The total 
one-sigma systematic error is shown for the oscillated prediction. The left panels show the MINOS-i- data 
on its own, and the right panels show the MINOS-i- data combined with the MINOS beam data. 


Amlj, 024 and 634 . Systematic uncertainties are reassessed, including beam systematic uncertain¬ 
ties which become more important due to the possible ND oscillations for large values of A/M 43 . 
Covariance matrices are used to include a total of 26 different systematic uncertainties in the fit, and 
the Feldman-Cousins approach is used to correct the log-likelihood surfaces [^]. For the results 
shown here, the full MINOS LE neutrino-dominated beam mode data set of 10.56 x lO^*’ POT is 
used. 

Fig. ^ shows a comparison of the FD data with the standard three-flavour prediction for the full 
MINOS LE neutrino-mode beam sample, including the corresponding Ear/Near ratios. 2712 v^- 
CC-like events, and 1221 NC-like events, are observed in the ED. After the final fit of the Ear/Near 
predictions to the Ear/Near data, MINOS gives the world’s strongest constraint on muon neutrino 
to sterile neutrino oscillations for <1 eV^. The corresponding contours can be seen in Eig. 
Moving forward, additional neutrino data from MINOS-h, and also data in antineutrino-enhanced 
beam mode will provide even better limits for the sterile neutiino disappearance analysis. 


4. Summary 

The MINOS experiment has been reborn as the new MINOS-i- experiment in the upgraded ME 
NuMI beam. It started taking data in September 2013 and has already contributed to the standard 
oscillations three-flavour analysis with additional atmospheric neutrinos data. This combined thi'ee- 
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Figure 6: Far Detector Spectra Comparison for the MINOS Sterile Analysis. This graph shows a comparison 
of the FD Charge Current and Neutral Current event data spectra in black to the MINOS three-flavour best-fit 
prediction in red. The red bands represent systematic uncertainties on the prediction. Various backgrounds 
are shown in colour. The Far/Near ratios for both data and prediction are also shown for both samples. The 
prediction is fit to the data Far/Near ratios to determine the oscillation parameters. 



8111 ^( 2024 ) 


Figure 7; Disappearance Exclusion for Sterile Neutrino Oscillations. The 90% C.L. exclusion limit in Ain^^ 
and sin^(2024) space is shown in black for the total MINOS neutrino exposure of 10.56x 10^'' POT. Other 
experiment results are shown in colour. MINOS sets the world’s best limit over several orders of magnitude 
in A77i43. 
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flavour analysis obtains the atmospheric parameter best-fit point of = 2 . 371 qq 7 x 10 ^ eV^ 
and sin^ 623 = 0 - 43 loo 5 in the inverted hierarchy phase space. The data shows a slight preference 
for the inverted hierarchy and the lower octant of 623 - A first look at the new MINOS-i- beam data 
is shown and is consistent with the updated results from the standard oscillations analysis. Finally, 
new MINOS results for the search for sterile neutrinos are presented and new constraints are placed 
on the possible phase-space accessible to potential sterile neutrinos using neutrino disappearance. 
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